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Positron annihilation spectroscopy is often used to analyze the local electronic structure of materi¬ 
als of technological interest. Reliable theoretical tools are crucial to interpret the measured spectra. 

Here, we propose a parameter-free gradient correction scheme for a local-density approximation 
obtained from high quality quantum Monte Carlo data. The results of our calculations compare 
favorably with positron affinity and lifetime measurements opening new avenues for highly precise 
and advanced positron characterization of materials. 
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The positron upon annihilation with its anti-particle, 
the electron, yields unique information about the elec¬ 
tronic structure of bulk materials 0,0 and nanostruc¬ 
tures 0 . The electron-positron density functional theory 
(DFT) 0 is used in order to obtain precise knowledge 
of the positron wave function and its overlap with the 
electron orbitals. The powerful combination of positron 
annihilation spectroscopy (PAS) and DFT calculations 
provides a highly accurate method for advanced charac¬ 
terization of materials [sj. 

Within the DFT framework, the generalized gradi¬ 
ent approximation (GGA) method to describe electron- 
positron correlation effects in solids has shown a system¬ 
atic improvement over the local density approximation 
(LDAlfor positron affinities and annihilation character¬ 
istics [6l-l9|. Until now, a dimensional analysis has been 
used to determine the form of the lowest-order gradient 
correction with a semi-empirical coefficient a. So far, the 
pragmatic approach has been to fit a to large databases of 
jositron lifetimes. Recently, both the LDA and the GGA 
IS El rave been improved on the basis of new quan¬ 
tum Monte Carlo (QMC) data for the electron-positron 
correlation problem in a homogenous electron gas (HEG) 
0 . 

However, one could claim that such good fits may be 
in some cases accidental [l^. Moreover, the present gra¬ 
dient corrections may also lead to some unphysical ef¬ 
fects in the electron-positron correlation potential near 
the nuclei; namely, its too large oscillations due to the 
shell structure of core electrons. Therefore, here we pro¬ 
pose to improve the GGA by extracting and deducing the 
a parameter from more fundamental physical principles. 
This more reliable derivation of a also reveals a gentle de¬ 
pendence of the local density reducing the gradient cor¬ 
rection near the nuclei. Thus, a becomes a function of 
the local density as well. 

In the case of the positron immersed in an electron gas 
the Goulomb attraction produces a cusp in the electron 


density at the positron site. The correlation potential 
describing the positron perturbation represents the elec¬ 
tronic polarization due to the positron screening and can 
be obtained via the Hellmann-Feynman theorem using 
coupling-constant integration as follows 0 


Vc{r) 




p{R) [g[r,R,Z)-l] 
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( 1 ) 


where p{R) [g{r, R, 2') — !] is the screening cloud density 
around a positive particle with charge Z {g{r, R, Z) is 
the particle-electron pair distribution function). The ef¬ 
fect of the density gradient on the correlation energy 
can be deduced from the distortion of the polarization 
cloud due to this gradient. For this purpose^one can 
use the dynamical structure factor ^(q,^) [lallgl of the 
HEG to show that in the high density limit the lowest 
order gradient correction is proportional to the param¬ 
eter e = (iVlnpl/qTp)^ (which depends on the ratio of 
the Thomas-Fermi length Atf = 1 /qtf and the inhomo¬ 
geneity length 1/|V lnp|). This correction is given by the 
expression 


^Vc{r) = P 
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( 2 ) 


where the constant /3 = 0.066725 Hartree is linked to the 
coefficient of the term in the density response func¬ 
tion wave vector expansion. The coefficient /3 has been 
calculated by Ma and Brueckner 0 and has been used 
by various authors 18 -2^. Eq. (I2|) is in fact similar to 
that used to compute the correlation en ergy [20j | for an 
electron gas with slowly varying density [21|. 

In order to interpolate to the case of rapid density vari¬ 
ations (i.e. large e), we use the formula 


Vc = exp(-ae/3) , (3) 

from Ref. 0 [see Eq. (7) there]. This formula is based 
on the scaling relation for the correlation potential, as 
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derived by Nieminen and Hodges (^ . But a is now 
a function of the local density (and thereby position). 
When we identify the first order expansion in e with the 
Ma and Brueckner’s result shown above, we find that 


to obtain calculated positron lifetimes within a precision 
of 0.1 ps and positron affinities within 0.01 eV. Here we 
consider only systems in which the positron density is 
approaching zero in the limit of an infinite crystal. 
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The quantity a remains a gentle function of the density in 
the valence electron region and at low density it becomes 
very close to 0. 05 12311 -- a value found earlier within the 
empirical GGA [lOj, [ll| . Interestingly, a happens also to 
be of the same order as the fraction Zc of an electron 
displaced in electron-electron correlation effects which is 
typically of the order of 1/20 of the electron charge 24, 


25|. 


Like the potential 14, the positron annihilation rate 
depends on electron-positron correlation effects and must 
be enhanced over the independent particle model. The 
electron-positron enhancement theory has some fea¬ 
tures in common with the interaction between a core hole 
and the conduction electrons treated both in X-ray emis¬ 
sion and in resonant inelastic X-ray scattering (^ . 
We can relate the correlation energy to the annihilation 
rate by using the scaling relation [2^ . Therefore, one ob¬ 
tains an electron-positron enhancement annihilation fac¬ 
tor 7 given by 


7 - 1 = ( 7 LDA _ gxp(-ae) , (5) 


The enhancement term 7 is used to calculate the to¬ 
tal positron annihilation rate or the inverse lifetime l/r, 
which is expressed through the simple relation Q 

i=7rrgc J d^r-/{r)p{r)\i}+{r)\^ , (6) 

where tq is the classical radius of the electron, c is the 
speed of light and is the ground state positron 

wave function. 

In this work, we have used the same accurate com¬ 
putational method described in Refs. [3 El- Elec¬ 
tronic structure calculations for selected materials were 
carried out using the self-consistent WIEN2k code [ 2 ^, 
which imposes no shape restrictions for the electron den¬ 
sity and the potential, while the positron wave function 
and energy were obtained using a Schrbdinger equation 
solver based on a finite difference method. The exchange- 
correlation potential for the electrons contains gradi¬ 
ent corrections within the scheme proposed by Perdew, 
Burke and Ernzerhof except in the case of the Ad and 
5d elemental metals since some of their calculated prop¬ 
erties (e.g. the lattice constant) become inappropriate 
when gradient corrections are used [s^. GGA correc¬ 
tions introduce cusps in the electron potential, negli gibl e 
in the LDA, which reflect the atomic shell structure [30j . 
Numerical parameters of the WIEN2k code as well as of 
the positron solver were tested and optimized in order 


TABLE I. Positron affinities (in eV) calculated according to 
various approaches; GC = original gradient correction with 
the Arponen and Pajanne potential [3i[ (q = 0.22), DB = 
Drummond et al. |l^ . DG = gradient correction with DB 
{a = 0.05), PF = parameter-free gradient correction with DB 
(varying a). The last column gives experimental values taken 
from Refs. The exceptions are C and Si (see Ref. 

0 and references therein) and MgO (Ref. [H). In the case 
of MgO an upper limit is given (see the text). 


System 

Structure 

GC DB 

DG 

PF Exp. 

Li 

Elements 

bcc -7.31 -7.02 

-6.95 

-6.96 

C 

diamond 

-1.33 -2.40 

-1.87 

-1.93 -2.0 

Na 

bcc 

-7.18 -6.89 

-6.80 

-6.81 

Al 

fee 

-4.21 -4.04 

-4.00 

-4.01 -4.1 

Si 

diamond 

-6.29 -6.47 

-6.33 

-6.35 -6.2 

Fe 

bcc 

-3.40 -3.76 

-3.62 

-3.67 -3.3 

Cu 

fee 

-3.76 -4.23 

-4.05 

-4.11 -4.3 

Nb 

bcc 

-3.61 -3.75 

-3.65 

-3.68 -3.8 

Ce 

fee, a-Ce 

-4.11 -4.16 

-4.07 

-4.09 

Ce 

fee, 7 -Ce 

-5.34 -5.31 

-5.23 

-5.25 

W 

bcc 

-1.72 -1.91 

-1.82 

-1.85 -1.9 

Pt 

fee 

-3.31 -3.77 

-3.61 

-3.67 -3.8 

MgO 

Compounds 

rock salt —5.56 —6.46 

-6.17 

-6.25 -5.2 

CU 2 O 

cuprite 

-5.88 -6.42 

- 6.21 

-6.26 

CeOa 

fluorite 

-6.55 -7.40 

-7.12 

-7.18 

YBaaCugOe 

tetragonal 

-6.11 -6.65 

-6.42 

-6.46 

YBaaCugOr 

orthorhombic 

-6.02 -6.78 

-6.52 

-6.58 

PrBaaCuaOr orthorhombic 

-5.81 -6.57 

-6.30 

-6.36 


DFT provides an excellent description of the Si elec¬ 
tronic structure both in the solid and liquid phases [38| . 
It is therefore natural to start our tests of the parameter- 
free GGA positron potential in Si. A meaningful observ¬ 
able to check is the positron affinity A defined as the sum 
of the electron and positron chemical potentials. In the 
case of a semiconductor, the electron chemical potential 
is taken from the position of the top of the valence band. 
Recently, Cassidy et al. [0| have shown that the tem¬ 
perature invariant time of flight (TOE) component for Ps 
emitted from the surface of p-doped Si(lOO) has a kinetic 
energy equal to 0.6 eV. This TOE feature is explained by 
a bulk positron picking up a valence band electron just 
beneath the surface to form Ps with a kinetic energy of 
K = Eps-I-A = 0.6 eV. Therefore, the experimental affin¬ 
ity for Si can be deduced to be A = —6.2 eV. When we 
use the GGA for both the electron and positron poten¬ 
tials, we find a theoretical value A = —6.35 eV, which is 
in excellent agreement with the value measured by Cas¬ 
sidy et al. while the corresponding LDA value shows a 
clear tendency to overestimate the magnitude of A. This 
LDA problem can be traced back to the screening effects. 
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TABLE 11. Positron lifetimes (in ps) calculated according to 
various approaches explained in the caption of Table m The 
last column gives experimental values discussed in Refs. [l^ . 
[TT| . The last experimental values for cuprates are extracted 
from Refs. (3^. [^ . 


System 

Structure 

GC 

DB 

DG 

PE Exp. 


Elements 





Li 

bcc 

283.2 

303.8 

316.2 

313.5 

291 

C 

diamond 

102.8 

94.6 

98.9 

97.7 

98-b 

Na 

bcc 

337.7 

343.0 

364.4 

360.5 

338 

Al 

fee 

154.2 

161.0 

164.1 

163.0 

160-t 

Si 

diamond 

222.7 

208.1 

217.3 

215.9 

216-t 

Fe 

bcc 

109.6 

102.1 

106.5 

104.7 

105-b 

Cu 

fee 

120.0 

107.4 

113.3 

110.9 

110 -t 

Nb 

bcc 

123.4 

120.9 

124.3 

123.1 

120 -t 

Ce 

fee, Q-Ce 

169.5 

165.0 

170.5 

169.0 

233 

Ce 

fee, 7 -Ce 

196.8 

194.1 

200.6 

198.9 

235 

W 

bcc 

102.7 

100.6 

103.4 

102.3 

105 

Pt 

fee 

105.2 

97.4 

101.3 

99.8 

99-t 


Compounds 




MgO 

rock salt 

146.2 

119.0 

128.5 

125.4 

130 

CU 2 O 

cuprite 

177.4 

147.3 

158.4 

154.8 

-174 

CeOz 

fluorite 

173.7 

138.2 

149.1 

146.0 

<187 

YBazCusOe 

tetragonal 

224.5 

175.4 

190.8 

186.5 

-190 

YBaaCusOr 

orthorhombic 

179.2 

142.4 

154.0 

150.5 

-165 

PrBa2Cu307 

orthorhombic 

180.4 

143.4 

155.0 

151.6 

-165 


In the GGA, the value of A agrees with the experiment 
by reducing the screening charge. Calculated positron 
affinities within LDA and GGA against the correspond¬ 
ing experimental values for different materials are shown 
in Table HI The trends follow those of Si, nevertheless 
the experimental values of A are often of earlier date 
and not always reliable. The corresponding positron life¬ 
times are presented in Table [III Clearly, the trends of 
the parameter-free GGA are very similar to the empiri¬ 
cal GGA [ISEl- In particular, one of the best result is 
given by A1 which was problematic in the original GGA 
scheme Positron lifetime measurements in Li and Na 
were performed before the advent of reliable spectrome¬ 
ters and fitting procedures, as discussed in detail in Ref. 
0 , and may be affected by significant errors. 

However, in the present scheme the positron has a 
slightly larger overlap with the core electrons as illus¬ 
trated in Figs. [T] and m for Si and Cu, respectively. Some 
noticeable jumps of e shown in Fig. [1] (d) and Fig. [2] (d) 
result in unphysically large local changes in the empiri¬ 
cal GGA correlation potential depicted in Fig. [T] (c) and 
Fig. [2] (c). These problems are now cured by the varia¬ 
tion of the function a in space illustrated by Fig. [T] (b) 
and Fig. |2| (b). Interestingly, a given by Eq. o seems 
to vary almost like the Thomas-Fermi length Atf and 
becomes very small close to the nuclei. Therefore, the 
cusps in the parameter free GGA correlation potential 
become more damped because of the reduction of the 
screening length in the core region. This effect is further 


documented by exp(—ae/S) factor plots (Fig. |T] (d) and 
Fig. |2| (d)) which define the reduction of the correla¬ 
tion potential in the core region. The exp(—ae/S) factor 
anticorrelates with the e parameter; i.e. a large inhomo¬ 
geneity corresponds to a small exponential factor. The 
variation of a in the core region should also improve the 
description of high-momentum annihilation spectra ob¬ 
served in coincidence Doppler broadening spectroscopy 
SEl and in angular correlation measurements (42| . 


The positron annihilation lifetime (PAL) provides a 
way to detect very small amounts of vacancy-defects in 
crystalline materials. Since thermalized positrons are 
trapped by vacancies before annihilating with electrons, 
their lifetime increases with respect the bulk values given 
the low electron density at the vacancy. For this reason, 
PAL has been widely used to characterize doped semi¬ 
conducting samples of silicon and other technological rel¬ 
evant materials 0. As shown by Table [III the positron 
bulk lifetime of Si is very well described by the present 
theory. Therefore, deviations from the theoretical life¬ 
time indicate the presence of imperfections in the sam¬ 
ple. In a post-silicon electronics era, engineered doping of 
oxide electronics, which is similar to conventional doping 
in semiconductor technology, offers much greater func¬ 
tionality including electronic control of redox chemistry 
with applications to batteries, photovoltaics and cataly¬ 
sis. In particular, a well characterized material is MgO, 
which is a simple binary oxide with rock-salt structure. 
In MgO, a magnetic moment can arise from the unpaired 
2p electrons at an oxygen site surrounding a cation va¬ 
cancy with each nearest neighbor oxygen carrying a mag¬ 
netic moment [d^. This magnetic property can be fine 
tuned to optimize spintronics devices. Goncerning PAL 
studies, Tanaka et al. 4^ have shown that MgO lifetime 
is significantly affected by Ga doping, which results in the 
creation of Mg vacancies. However, when the number of 
Mg vacancies decreases the lifetime converges to the bulk 
value 130 ps [i^, which is in reasonable agreement with 
the present theory. A reliable experimental TOF study 
of MgO 3^ reports a Ps emission peak energy of 2.6 eV. 
Since Ps is already formed in the bulk of MgO, the kinetic 
energy is given in this case by K = Eps + A — Eb + Eg, 
where As is the Ps binding energy inside the MgO ma¬ 
trix and Eg = 7.8 eV is the energy gap of MgO. Using 
our calculated affinity, we deduce that Eb = 5.75 eV, 
which is consistent with typical values of Ps binding en¬ 
ergy in the bulk [i^. In fact, this value must be smaller 
than Eps because of screening effects in the bulk. 

Ceria T7| is another oxide which has attracted consid¬ 
erable interest because of its applications in solid oxide 
fuel cells. It can be noted that by removing all the oxy¬ 
gen atoms, one recovers the fee structure of Ce. Exper¬ 
imentally, positron seems only to detect the 7 phase of 
Ce because of its stronger affinity with respect to the a 
phase. Interestingly, the experimental ceria lifetime 189 
ps [ 4 ^ appears to be much closer to theoretical value of 
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FIG. 1. (Color online) One-dimensional profiles of (a) the 
electron density (including atomic orbitals), (b) the a param¬ 
eter (and the Thomas-Fermi length Atf); (c) the positron cor¬ 
relation potential, and (d) the exponential factor exp(—ae/S) 
(and e parameter) along the [100] direction in Si for LDA 
(DB), the empirical (DG) and the parameter-free (PF) GGA 
approaches. Si atoms are located at 0 and 10.26 au along 
[100 ]. Atf Ei-nd e are shown for the purpose of observing cor¬ 
relations with corresponding quantities (the scales of Axf ^nd 
e are different than those for a and exp(—Qe/3), respectively). 

7 -Ce rather than ceria. A possible reason for this dis¬ 
crepancy is that real samples can always contain patches 
of 7 -Ce which strongly attract the positron because of 
their higher positron affinity. In this context, we should 
keep in mind that oxygen is very mobile in ceria. 

As an example of advanced characterization, we now 
show that positron annihilation spectroscopy can be use¬ 
ful to understand the role of oxygen-related defects in 
high temperature superconductivity 14911. In practice, by 
comparing the experimental lifetimes |36l| to an accurate 
theory it is possible to deduce that positrons are trapped 
at oxygen vacancies in the superconducting compound 
YBa 2 Cu 307_5 while this trapping becomes negligible in 
the non-superconducting compound where Y has been 
replaced by Pr. When positrons become completely de¬ 
localized for temperatures higher than 400 K, the life¬ 
time becomes almost identical in the YBa 2 Cu 307 and 
PrBa 2 Cu 307 compounds in agreement with our calcu¬ 
lations reported in Table m Moreover, the calculated 
lifetime in the tetragonal YBa2Cu306 lattice is 36 ps 
longer than in the orthorhombic YBa^C^O?. Such dif¬ 
ference is consistent with experiments [ST]. Curiously, 
the calculated positron affinity seems to indicate that Ps 
is emitted with about 0.15 eV higher kinetic energy from 
YBa 2 Cu 306 and PrBa 2 Cu 307 than from YBa 2 Cu 307 . 


FIG. 2. (Golor online) One dimensional profiles for Cu as 
explained in the caption of Fig. [T] Gu atoms are located at 0 
and 6.83 au along the [100] direction. 


Nevertheless, since the present DPT calculations fail 
in describing the insulating phase of YBa 2 Cu 306 and 
PrBa 2 Cu 307 we should take the positron affinity calcu¬ 
lated values for these two compounds with caution. 

In conclusion, we have demonstrated that the 
parameter-free GGA truly provides a simple, yet accurate 
step beyond LDA. It is also reassuring that the most re¬ 
liable electron-positron LDA parametrization (based on 
the QMC simulations) combined with the parameter free 
gradient correction gives the best results compared with 
any of the older LDA potentials. Further studies com¬ 
bining the present approach with well-converged momen¬ 
tum densities calculations [H^ are needed to check if first 
principle methods can soon improve the agreement over 
empirical approaches (42| . 
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